Droplet spreading and pinning on heterogeneous substrates 



Mikko J. Alava 
Aalto University, Department of Applied Physics, 
PO. Box 14100, 00076 Aalto, Finland 

Martin Dube 

Centre de Recherche sur les Materiaux Lignocellulosiques, 
Universite du Quebec a Trois- Rivieres, Trois- Rivieres, Canada C9A 5H7 

Abstract 

The contact angle of a fluid droplet on an heterogeneous surface is analysed using the statistical 
dynamics of the spreading contact line. The statistical properties of the final droplet radius and 
contact angle are obtained through applications of depinning transitions of contact lines with non- 
local elasticity and features of pinning-depinning dynamics. Such properties not only depend on 
disorder strength and surface details, but also on the droplet volume and disorder correlation 
length. Deviations from Wenzel or Cassie/Baxter behaviour are particularly apparent in the case 
of small droplet volumes and small contact angles. 
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I. INTRODUCTION 



Liquid droplets on a solid surfaces pose several interesting theoretical and experimental 
challenges and has a direct importance in several industrial processes. In general, the inter- 
action of the liquid with the surface is characterized through the contact angle that a droplet, 
with known volume, makes with the surface. For a flat, chemically homogeneous surface, the 
equilibrium contact angle is easily obtained from an energy minimization procedure. The 
result, known as Young's law, includes the three surface energies (solid-liquid, solid-vapor, 
liquid-vapor, •jgi, 7^^, 7, respectively) with a "weight factor" for the last one, the cosine of 
the contact angle, cos^. The dynamical approach to equilibrium is also well understood, 
both for complete jll and partial ^ wetting. In reality, the static contact angle, as well as 
its dynamics, is obviously influenced by the variations in the physicochemical or geometrical 
properties of the substrate on which the droplet spreads l3|-l5|] . Such surfaces are present 
in an enormous variety of industrial processes. Typical examples range from printing and 
coating to painting, as well as the creation of tailored superhydrophobic surfaces [6|. 

The static contact angles for surfaces with local variation can be obtained from variants of 
the original energy minimisation. The classic Cassie-Baxter treatment for spatially varying 

n 

surface energies, considers a simple spatial average, ie., (7(r)) [7]. Likewise, Wenzel law for 
surfaces with height variations introduces the ratio - always larger than unity - of the total 
surface to the projected surface as a correction factor [8]. Clearly both Cassie-Baxter and 
Wenzel treatments may not apply in practice, since a droplet spreading towards equilibrium 
will encounter several pinning states in which it can be trapped, thereby never reaching the 
theoretical result [o], [lO] . This is a much-studied issue both experimentally and theoretically. 



particularly for engineered surfaces. 



Ill 4l4| . Approaches that have been tried range from 



hydrodynamic simulations of spreading to quasi-static considerations [l5|, [l6|. In all cases, 
a crucial concept is the spreading power, or the difference of the actual contact angle from 
the static or equilibrium one. 

Here we analyze the consequences of surface heterogeneity on the spreading and final 
static state of liquid droplets. Our approach centers on the evolution of the droplet radius, as 
defined by the contact line that separates the wetted and non-wetted parts of the substrate. 
The radii are then easily related to values of contact angle through volume conservation. 
The droplet radius evolves under the influence of three factors: i) the imbalance in surface 
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tension, directly related to the spreading power, ii) the quenched noise, induced by the 
locally varying surface properties, that introduces angular deformation in the radius, and 
iii) contact line elasticity that tends to smoothen radius angular variations. 

We use contact line pinning and depinning dynamics to follow the evolution of the droplet 
radius. In the initial stages of spreading, surface tension imbalance controls the dynamics. 
However, as the droplet spreads, disorder effects start to dominate and the contact line comes 
to a halt before an equilibrium state can be reached. Disorder, elasticity and spreading thus 
predict a pinning transition for the contact line, and thus a droplet pinning radius that is 
sample-dependent. We show how this depends on the the independent parameters of the 
theory: disorder strength, disorder length-scale, droplet radius and volume, and the fluid 
properties such as surface energy and contact angle. 

We further obtain the statistical distributions of the pinning radius. Small droplets are 
strongly effected by disorder, leading to results that are markedly different from what could 
be expected within a Wenzel or Cassie-Baxter framework. Disorder effects are apparent in 
the case of droplets with a small contact angle and disorder variations that are correlated 
on a short length-scale with respect to the droplet radius. The pinning predictions for the 
contact angle have simple experimental consequences: they can be tested by varying for a 
given type of surface the nominal contact angle (by changing liquid for instance) and by 
varying the droplet size. 

The rest of this paper is structured as follows. In Section Ull we outline the basic spreading 
phenomena as well as the approach to the pinning phase. Droplet spreading in the pinning 
phase is analyzed in details in Section and Section HV] finishes with conclusions. 

II. DROPLET SPREADING AND APPROACH TO PINNING 

We consider a liquid droplet with air/liquid interfacial tension 7, density p and viscosity 
r] deposited on a macroscopically fiat surface. When the size of the droplet is smaller than 
the capillary length {'j/pgY^'^ ~ 1 mm, the droplet has the overall shape of a spherical cap 
with basal radius R, the air/liquid interface joining the solid surface with contact angle 6. 
For a perfectly homogeneous surface, the equilibrium contact angle 6eq can be obtained from 
Young's law: 

7COs6'eq = 7s„ - 7si, (1) 
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related to the equilibrium radius Req via the drop volume Vt: 

n = ^ (2-3cose + cos^e) . (2) 

3 sin^ ' ^ ' 

This relation reduces to = ttR'^O/A for small contact angles 0^1. 

Equation [1] has to be supplemented by a correction term involving the line tension r, 



related to k, the curvature of the contact line 
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7" 

COS 6 = cos 6eq H K (3) 
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The actual value of the line tension remains a difficult quantity to measure, due to chemical or 
topographical heterogeneity of the substrate and the difficulty to perform experimental work 
with sub-micron sized droplets. In general, the value of the line tension can be estimated 
to be r ~ 10^^'^ A^, thus representing a relevant correction only for droplets with nanometer 
range dimensions. 

Substrate disorder induces variations in the contact line and hence in the local values of 
the contact angle and radius of the droplet. If x is the local coordinate along the contact 
line, the local droplet radius is R{t) + h{x,t), where R{t) is the spatially averaged radius. 
Variations along the contact line are connected to variations of the local contact angle 6{x, t), 
since the overall volume of the drop is conserved. To first order in the radius variations, the 
variations in the contact angle are Q 

eix,t)=9it)^l + ldx'^^^^ (4) 

Correction terms coming from the overall curvature of the droplet can be neglected for small 
variations of the contact line. 

Upon deposition on a surface, a droplet will tend to its equilibrium shape through hydro- 
dynamical spreading, under the influence of uncompensated Young's force Sy{9{x,t)) and 
disorder variations 5d(6'(x, t)) |l| 

where / is a numerical factor arising from finite slip at the contact line 

For small values of the contact angle 6 <^ 1 and to first order in contact line deviations 
(using Eq. H])), the uncompensated Young force, corrected with a term arising from the 
dissipation at the contact line, 

Sy{9{x, t)) = tan(^(a;, t)) (cos(^eg) - cos(^(a;, t)))) (6) 



is decomposed into an overall term independent of contact line variations and a non-local 
term 

Sy{e{x, t)) = s{e{t)) + s,i{e{t), h{x, t)) (7) 

where the spreading power 

s{e) = e{e'-e%), (8) 

is balanced by the elastic restoring force 

Sei{9, h{x, t)) = 9{39' - 6%) I dx'^^^^ (9) 

For simplicity, we consider a chemically disordered substrate, which contributes a 
quenched random force. Again, in the limit 6 <^ 1 and to first order in the contact line 
variations, 

Sd{e) = %{h{x,t)) (10) 
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where 

6'y{x, h{x)) = Sjsvix, h{x)) - 6'ysi{x, h{x)). (11) 

The noise is generally delta-correlated over two microscopic scales ^y, ^x, parallel and per- 
pendicular to the direction of spreading, dependent on the nature of the substrate. It 
includes the chemical contrast (via the variation of the surface energies, and 7^^, where 
the indices denote the solid-vapor and solid-liquid cases, respectively) and may also include 
surface roughness through local surface tilts [l|, [2^ . In what follows, we consider the slowly 
evolving dynamics of the drop. This can be achieved experimentally by slowly depositing a 
droplet on the substrate, thereby avoiding initial impaction effects. After that, the viscous 
contact line dynamics dictates the (slow) evolution of the droplet [2l| . The analysis does not 
apply to volatile liquids and the Ohnesorge number Oh = ri/{pR'y) is an appropriate param- 
eter to estimate how viscous the liquid must be in order to prevent inertial as well as front 
instability (coming from evaporation) effects to occur 22]. We only consider microscopic 
disorder such that the length scales and are the smallest length scales of the problem. 
The droplet then globally keeps a spherical shape, with disorder inducing variations in the 
contact line itself. 

On a disordered surface, the equilibrium contact angle of a droplet may be markedly 
different from the prediction of Eq. ([T]), depending on the type or disorder. For example, 
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bi-modal chemical disorder, where a fraction / of the surface has surface tension 7 + A7, is 
usually analyzed in terms of Cassie-Baxter equation 

.A7 



cos 9 



[CB] 



COS 9eg - f 
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(12) 



which is essentially a weighted average of the equilibrium contact angle on each surface. On 
the other hand, topographical disorder (surface roughness) is encompassed within Wenzel's 
result 



cos 9 



w 



Ar COS 9, 



eq 



(13) 



where the relative area Ar > 1 is the ratio of the total surface to the projected surface under 



the drop (see 23|,|2J] for recent similar analyses). It is interesting to note that, through Wen- 



zel's analysis, roughness increases the wetting or non-wetting tendencies originally present 
in the problem, ie., if 9eq < tt/2, 9y/ < 9f.q until cos9eq = at which point a wetting film 
should be formed. 

Both approaches however neglect that spreading is a dynamical phenomenon and that a 
droplet may become pinned in several configurations before reaching equilibrium. Spatial 
scales and correlations of the disorder are also not included. 

During spreading, the contact line roughens due to the quenched disorder of the substrate, 
a phenomenon analyzed extensively in the literature 
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25l-l27l|. also in the context of the 



hysteresis for advancing and receding contact lines 28|]. When disorder becomes relevant the 
interface dynamics enters a critical regime where pinning occurs, and the interface propagates 
through avalanches. This regime is characterized by a force Fc, at which an interface becomes 
pinned. Around F^, the contact line develops local self-affine fluctuations (roughness) w ^ 
where r is the length of the interface and ( ~ 0.38 the roughness exponent specific to contact 



line motion 
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H3l|. 



The propagation of the interface by avalanche is characterized by a series a correlation 
lengths and critical exponents. During the avalanche, a portion of the avalanche, with lateral 
size ^, moves by a distance w ~ The duration r ~ of the interface is characterized by 
the dynamical exponent z. The lateral extent of the avalanches is a correlation length related 
to the driving force F and the pinning force Fc by the critical exponent u: C, ^ {F — -F^)"^. 
Finally, the velocity of the interface scales as f ~ (F — F^)^. Critical scaling implies that 



(3 = z/(z — () and, for non-local e 
of contact line motion (3 ~ 0.625 



astic interface, z/ 



() ^. Again, for the specific case 



For any finite system, and thus the corresponding critical angle 9c have a finite-size 
correction and in particular a sample (and disorder) -dependent actual critical value, with 
an universal probability distribution. This distribution is characterized by its width which 
decays with the interfacial length and depends only on the stren gth of the disorder, measured 
by the prefactor of the disorder two-point correlation function 33 1. 

Although a complete solution of roughening requires additional non-linear terms, Eq. 



([5]), together with typical scaling arguments 
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26| . already provides much information. 



Balancing the elastic restoring force (Eq. (|9])) against disorder (Eq. (fTOj) ) yields a length 
scale Ic = ^±(36*^ — O^^Y'^/T where T = (57)^/7 describes the pinning strength of disorder. 
For length scales / ^ Ic, the elastic restoring force dominates while disorder dominates for 
/ ^ Ic- The contact angle 6p at which pinning first becomes relevant is then obtained by 
balancing the spreading power (Eq. ([8])) against either the elastic or pinning force at the 
length scale Ic- To first order in {9p — 9eq)/9eq, 

Or> = Oeq + ^^ (14) 
Through volume conservation (using Q = 7tR^9/A), this translates to a pinning radius 

J^ = _L + ;™,i,3^£Z.. (15) 

The volume of the droplet thus plays a crucial role. Large droplets will enter the pinning 
regime already close to equilibrium while, for smaller droplets, the ratio Rp/ Reg can be much 
smaller than 1, even more so for strong wetting 9eq '^l- 



III. SPREADING IN THE PINNED PHASE 



We now consider spreading on a disordered substrate. The case of a contracting droplet 
(for equilibrium contact angles in excess of 7r/2) is an easy extension not reported here. The 
droplet has initial basal radius Rq <^ Reg and increases its radius through hydrodynamical 
spreading until Rp- Then, disorder becomes relevant and spreading continues if the spreading 
power (cf., Eq. ([8])) at position R, S{R > Rp), can overcome the combined pinning potential 
coming from the disordered substrate and the elastic restoring force, denoted Sc- At the 
first pinning radius, the spreading power 



where the last form is in the hmit {Req — Rp)/Req ^ 1 with 5*0 = QOl^/Req- Sp corresponds 
to a contact angle value 6p. In the similar limit, the elastic and disorder- induced forces can 
be written 



dx' 



I h{x',t) 



+ 



eq 



(17) 



""7 (x - x'y 

An analysis of spreading to a radius past Rp necessitates the knowledge of the pinning 
force distribution resulting from Eq. (^^. Extensive numerical simulations of a driven 
interface subjected to the combined elastic/disorder force (Eq. f|T7j) ) have shown that the 
probability F{S{r)) of propagating past a radius r derives from a distribution of pinning 
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forces for an interface of length r, f(sc,r) 

F{S{r))= rdsj{s,,r 
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f{sc,r)dsc 



The distribution follows a scaling form 



1-C 



Sn 



Sn 



1-C 



with a scaling function ip{f) independent on the details of the disorder 



(18) 



(19) 



Close to 

pinning, {Sp — Sc)/Sp <^ 1, the scaling function has a power law behavior, ip{f) ~ /'*' with 

7 = C/(i-C). 

In this critical regime, the likelihood P{R) that the droplet radius will at least be R 
is obtained from a probabilistic argument, an approach also used for indentation crack 
propagation [34]. The motion of the droplet (see Fig. [1]) consists of a succession of steps 
over independent configurations of the pinning disorder - ie. the distance covered is divided 
into uncorrelated increments. The size of the steps can be inferred from the critical dynamics 
of the contact line motion. After step i, during the sequence of jumps, the pinned contact 
line explores the combined potential (elastic plus disorder) over a distance w{R), finding a 
configuration that tend to minimize its energy. Upon depinning (ie., step i + l), the contact 
line moves, by a distance w{R) 25|, |26|, into a new configuration. At this point, both the 
contact line configuration as well as the combined potential are completely uncorrelated with 
the previous one. The appropriate step size during the sequence of events is thus w{R), the 
roughness of the contact line. A droplet will thus reach a given R{t) if it has passed through 
all the previous pinning zones without being stopped. The probability of such a chain of 
events is 

P{R) = n FiSiRp + z/A)) > S,) (20) 
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where F{S{r) > Sc) is the probabihty that the spreading power at radius r is larger than 
the pinning spreading power Sc- Time is not explicitly included in this argument. In the 
continuous limit, Eq. (12U]) can be rewritten 

-R 



P{R) ~ exp 



log(F(5(r))) A(r)dr 



(21) 



where the zone size 1/A is related to the droplet radius and the correlation lengths of the 
disorder, ^ oc 

For a droplet size large as compared to the scale of heterogeneities, effective contact angle 

values will always remain in the immediate vicinity of 6p and the power-law behavior of 

1 

at the origin can be used to approximate log(F(S'(r))) oc — ( ^5~^ ) ^ which leads to 



where the prefactor A oc . depends only on the material parameters through the corre- 
lation lengths of the disorder C,\\ and C,j_- The use of the linearized form for the spreading 
power, Eq. (fT6|) and a change of variable x = R/Rp finally yields: 

P{xRp) ^ exp ^ duu^''^{u - 1)1/(1-^) I , (23) 

with B = Req/Rp — 1. The integral in Eq. (I23ll exhibits an universal form which only depends 
on A, a disorder scale parameter, and Rp, related to Reg and the strength of the disorder 
through Eq. f fTSj) . The size of the droplet is thus implicitly present in Eq. f l23|) . 

To analyze the results, it is convenient to set C,± = = C,- The relevant dimensionless 
ratios are then Rp/R^q, a measure of the influence of disorder strength, droplet volume 
and equilibrium wetting properties, and Req/^, which relates the droplet typical size to the 
spatial structure of the disorder. These ratios are made apparent from the limit x ^ 1 of 
Eq. f l23|) which reads 

The probability for the droplet to reach a given R thus decays quickly close to -Rp. This 
decay is sharper for drops that are large compared to the disorder scale (-Rp ^ ^) than 
for smaller drops. This tendency to cluster around -Rp drastically increases as the ratio 
Req/Rp — > 1, which occurs for weak disorder or very large drops. 



The importance of the ratio Rp/^ is clearly shown in Figure 121 where the probability to 
reach a radius R as calculated from Eq. fl23l) is shown for various ratios Rp/^ and Rp/Req. 
For values Rp/^ ^ 1, this probability drops sharply and the drop remains essentially pinned 
at a radius Rp. It is only for relatively small values of this ratio that the probability to 
reach a radius larger than Rp increases significantly. In other words, occasionally for small 
droplets the spreading can get closer to the equilibrium radius. 

The effect of disorder on the values of contact angles at pinning can then be elaborated 
by comparing the predicted value to 9eq using volume conservation, Eq. (|2]). Figure [3] shows 
the average contact angles calculated from the theory for two values of 9eq- At large values 
of the ratio Req/C,, the final value of the contact angle is essentially determined by the value 
of Rp, ie., strong disorder, characterized by a small ratio Rp/Req leads to a larger apparent 
contact angle. In such a case, it is only for a small value of the ratio Req/C, that the apparent 
contact angle can be quantitatively close to the expected 9eq. 

In Figure m we show for three different disorder (F) values how Rp/Req itself scales for 
various contact angles and a fixed ^ taken to be 10 micrometers. The values of F are taken 
to have representative values; note that the parameter measures the relative variation of 
the surface energy due to chemical disorder or due to roughness. Out of a variety of cases, 
we depict nine representative ones to show the trends. Figure shows the actual droplet 
volumes at hand for each contact angle and disorder, while Figure Hb shows the ratios of the 
pinning radius to the equilibrium one. The trends themselves are obvious (larger contact 
angles lessen the effect of pinning, while stronger disorder works to the other direction), and 
the final prediction is then to be computed similarly to Figure [3] for all the cases, separately. 

It is thus interesting to note that the wetting properties of the surface, present through 
6eq also enter explicitly the problem through the values of Rp and Req. This is in contrast 
with Wenzel law or Cassie law which relate the apparent contact angles to the equilibrium 
contact angle through a set of constants independent of the nature of the surface or of the 
liquid. For naturally wetting surfaces {cos 9eq > 0), Wenzel law predicts that the apparent 
contact angle is larger than the equilibrium contact angle. This however neglects the fact 
that spreading is a dynamical process and that pinning of the contact line impedes the 
droplet from reaching an equilibrium state. 

The temporal aspects of the radial pinning process can be obtained from the relation 
between the interface velocity and the driving force in the critical regime. Hydrodynami- 
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cal spreading of the droplet occurs until the droplet reaches the radius Rp, at time Tp ~ 
{Rp/Reqf{Rp'r]/'^)/9l . At this point, the spreading velocity Vp = V{Rp) = '^9p{9'i — 9eq)/r7]l. 



After this point, motion proceeds by avalanches, with a velocity 



V{R) = V{Rp 



Req — R 



Req — R 



(25) 



where again /3 ~ 0.62. The time needed to move across a given shell of thickness 
A"^(r) is simply AT(r) = X^^{r) /V{r) and the total time to pinning is obtained from 
summing the successive contributions of each shells. In the continuous limit, T{R) = 



Tp + 



[Req-Rf CR 



f 

JR 



iRp {R^^Pi^ 

T{R) = Tp + - 



dr which is calculated to 



1 Req 


( Req] 
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(26) 



and can be reduced to T{R)—Tp ~ {R — Rp) /Vp in the limit of strong pinning {Rp/R^q ^ 1). 
This result, averaged over the distribution of pinning radius P{R) is illustrated in Fig. [5] 
and shows that that the time to reach a given R > Rp increases as a function of Rp/Req, 
particularly for R close to Reg. 



IV. CONCLUSIONS 



To summarize, we have presented a statistical physics theory of droplet asymptotic con- 
tact angles on heterogeneous surfaces. This allows us to identify the important quantities, 
such as the first pinning radius Rp and the ratios to the disorder scale and the equilibrium 
radius, Rp/R^q and Rp/^. The theory presented here applies to droplets with typical lengths 
above the nanometer scale, such that line tension effects are irrelevant. It is also important 
that radius variations are smaller than the radius of the droplet, for it to keep a compact 
circular shape on average. This last requirement implies that the disorder length scale 
^ < Rf.q- The contact line dynamics model loses its validity in certain cases - one example 
is when long-range microscopic forces make it invalid, or when the presence of features such 
as corrugations make it so that the coarse-grained surface tension indirectly assumed is not 
present. 

Our results predict a dependence of the average contact angle on the volume of the 
droplet. At constant disorder strength, we expect small droplets to exhibit a markedly 
larger contact angle at pinning than larger droplets, even more so when the droplet is large 
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with respect to the spatial scale of the disorder. The predictions of the theory can be easily 
tested through repeated spreading experiments using droplets of different volumes on the 
same substrate, so that the first pinning radius Rp (cf., Eq. fllSp ) only depends on the drop 
volume. It is also possible to test the theory using different liquids on a given surface. Large 
variations of the apparent contact angle due to the final stage of spreading are expected for 
small droplets. Many of the consequences of the theory of elastic manifolds are in contrast 
to static, energy-minimization based results as the Wenzel or Cassie-Baxter laws. 

Further theoretical developments include extending the theory presented in this paper to 
receding radii on hydrophobic surfaces and to develop similar arguments for a structured 
surface on which gas phase pockets can develop. It is also clear that the theory can be 
tested through large scale numerical simulations of the relevant hydrodynamics equations. 
Extending the probabilistic argument to the finite-temperature very-long time creep motion 
regime, relevant for contact lines as elastic manifolds, is also possible. 
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-program) and TEKES FinNano program via the Silsurf -project. MA would like to thank 
for the hospitality of prof. E. Frey and the Arnold Sommerfeld Center at the Ludwig- 
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FIG. 1: Schematic of the motion through jumps over independent pinning shehs of a size (rough- 
ness) u; ~ i?'-. At each shell, the local critical force follows the distribution of Eq. ()19p . The global 
driving force decreases as the spreading proceeds. 
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x = R/R 

eq 

FIG. 2: Probability for the droplet to reach a given radius R. The probability P{R) = 1 until 
R = Rp after which it decays depending on the values of the ratios Rp/Req and R(,q/^. Solid, 
dotted and dashed lines respectively correspond to Req/^ = 25, 100 and 400, while two values, 
Rp/Req = 0.8 (left set of curves) and Rp/Req = 0.9 (on the right) are used. 



15 



110 



o 

^ 80 

70 



T 1 1 — I — TT" 



T 1 1 1 1 — I — TT" 







O 


0.8 


100 




o 


0.9 






□ 


0.95 


90 









o 



o 



o 



o o 



o o 



1 1 1 — I — I — I — n~| — 

o o o o o o 



a) oo^nnaa°°°°°°°° 
^ n n u 



□ □ 

J I I I I I I 



J I I I I I L_L_ 



_l I I I I I L_L 



10 



100 

R 

eq ^ 



1000 



55 



a 

o 



25 



T 1 1 — I — TT" 



T 1 1 1 1 — I — TT" 



T 1 1 — I — I — I — n~ 



50 




• 0.8 




♦ 0.9 


45 




■ 0.95 


40 




• 


35 




• 


30 


^ b) 





♦ ♦ 



♦ ♦ ♦ ♦ 



♦ ♦♦♦♦♦ 



J I I I L_L_ 



_l I I I I I L_L_ 



J I I I L_ 



10 



100 

R 

eq ^ 



1000 



FIG. 3: Average values of the contact angle at pinning for a) ^gg = 75° and b) 6f,q = 30°. The 
contact angle is shown as a function of ReqH for three ratios Rp/Req = 0.8, 0.9 and 0.95. 
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FIG. 4: Volumes of droplets in nanoliters as a function of contact angle (in degrees, the values are 
from 5.7° to 18.1°) and disorder strength T (a)), b): the resulting values of Rp/Req. The disorder 
scale ^ is set to 10 fim 
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